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1 Introduction 

Trapped Bose-condensed atomic gases Si are remarkable because, in spite 
of the fact that these are very dilute systems, they exhibit robust coherent 
dynamic behaviour when perturbed. These quantum "wisps of matter" are a 
new phase of highly coherent matter. While binary collisions are very infre- 
quent, the coherent mean field associated with the Bose condensate ensures 
that interactions play a crucial role in determining the collective response of 
these superfluid gases. 

In our discussion of the theory of collective oscillations of atomic con- 
densates, which is the main theme of these lectures, the macroscopic Bose 
wavefunction <£>(r,i) will play a central role. This wavefunction is the BEC 
order parameter. The initial attempts at defining this order parameter began 
with the pioneering work of Fritz London □ in 1938, was further developed 
by Bogoliubov □ in 1947 and finally formalized in the general quantum field 
theoretic formalism of Bcliaev U in 1957. The first extension of these ideas 
to inhomogeneous Bose condensates was by Pitaevskii and, independently, by 
Gross in 1961, which led to the now famous Gross-Pitaevskii (GP) equation 
of motion for $(r, t). Most of this early work was limited to T = where, in 
a dilute Bose gas, one can assume all of the atoms are in the condensate. 

In Section |2|, I will first review the dynamics of a pure condensate at 
T = 0, based on solving the linearized GP equation. This limit is especially 
appealing since one can ignore all the complications which arise from the 
presence of non-condensate atoms (the thermal cloud). We will discuss the 
normal mode solutions of the T — GP equation of motion using the "quan- 
tum hydrodynamic" formalism, which works in terms of the local condensate 
density n c (r,i) and supecfluid velocity v c (r,t). Within the Thomas- Fermi 
approximation, Stringari □ has shown that the equations of motion for these 
two variables can be combined to give a wave equation for oscillations of the 
condensate. 
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In Section [|, we derive a generalized form of the GP equation for $(r,i) 
which is valid at finite temperatures. 13 It involves terms which are coupled 
to the non-condensate component (the thermal cloud) and thus its solution 
in general depends on knowing the equations of motion for the dynamics 
of the non-condensate. In the present lectures, we will restrict ourselves to 
finite temperatures where the non-condensate can be described by a quantum 
kinetic equation for the single-particle distribution function /(p,r,i). In this 
Boltzmann equation, the relatively high energy non-condensate atoms are 
simply free atoms moving in a self-consistent Hartree-Fock mean field. 

A unique feature of a Bose-condensed gas is that the kinetic equation 
for /(p,r,i) involves a collision integral (denoted by C\% [/,$]) describing 
collisions between condensate and non-condensate atoms. The generalized 
GP equation for $(r, t) also has a term which is related to the Cn collisions. 
This gives rise to damping of condensate fluctuations. In Section |[ we derive 
a finite T Stringari wave equation 13 with damping using the static Popov 
approximation. This means the thermal cloud is treated as always being 
in static thermal equilibrium, with /(p,r,t) = /o(p ; r ) being given by the 
equilibrium Bosc distribution. 

In Section we turn to a detailed treatment of the dynamics of the 
coupled condensate and non-condensate components, starting from the finite 
T generalized GP equation for Q(r,t) and our kinetic equation for /(p,r,i). 
We derive a new set of two-fluid hydrodynamic equations, following the recent 
work by Zaremba, Nikuni and Griffin. 13 This derivation assumes the non- 
condensate is in local hydrodynamic equilibrium, induced by rapid collisions 
between the atoms in the thermal cloud. As a result, the non-condensate is 
completely described in terms of the non-condensate density n(r, i), velocity 
v„(r,i) and local pressure P(r,t). Both the condensate and non-condensate 
exhibit coupled coherent oscillations at the same frequency. 

In Section ||, we review the famous two-fluid hydrodynamic equations first 
derived by Landau in 1941 and which form the_hasis of our understanding of 
the hydrodynamic behaviour of superfluid 4 He. l§[L2l We discuss first sound and 
second sound and consider how they differ in superfluid 4 He and in a uniform 
Bose gas. We also point out the appearance of a new hydrodynamic zero fre- 
quency mode which appears in our generalized two-fluid equations. This mode 
is associated with the relaxation time for the condensate and non-condensate 
atoms to come into "diffusive" equilibrium with each other. Inclusion of damp- 
ing due to thermal conductivity of the thermal cloud shows that this mode is 
the analogue of the well-known thermal diffusion mode in a classical gas. 

We close this Introduction with some general references which may he 
useful to the reader. The 1998 Varenna Summer School lectures on BEC □ 
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has excellent articles on recent research on BEC in atomic gases, including 
very detailed reviews of the experimental work carried out by the JILA 1121 
and MIT lij groups. For a general introduction to the theory of trapped Bose 
gases, we recommend the Reviews of Modern Physics EJ article by a lead- 
ing theoretical group at the University of Trento. This authoritative review 
concentrates on thermodynamic properties and the collisionless dynamics. A 
recent long paper by Zaremba, Nikuni and the author □ discusses the details 
of the derivation of the coupled hydrodynamic equations for the condensate 
and the thermal cloud. The present lectures may be viewed as an introduction 
to this paper, with more emphasis on the general structure of the derivation 
and the implications of this new set of superfluid equations. 



2 Pure condensate dynamics (at T = 0) 

This section will be largely a review of standard material but will provide 
the starting point for generalizations in the following sections. As we have 
mentioned, at T — 0, one can assume that all atoms in a dilute Bose gas 
are described by the macroscopic wavefunction $(r, t). This obeys the time- 
dependent Hartree equation of motion first written down 40 years ago by 
Pitaevskii and Gross, 



<9$ 



2m 



+ V ex (r) + V H (r, t) 



Here the trap harmonic potential is 



,(r) = imcjgr 2 



(isotropic) 



and the self-consistent condensate Hartree potential is 



*(r,t). 



(1) 



(2) 



V H (r, t)= J dr'v{r - r> c (r, t) = gn c (r, t). (3) 

As usual, since we are interested in extremely low energy atoms, we use 
the s— wave approximation and approximate the interatomic potential by the 
pseudopotential 

. Aitah 2 , ,. 

v(r — r j = d(r — r ) 

m 

- gS(r - r'), (4) 

where a is the correct s— wave scattering length. For further discussion of the 
effective interaction to use in the Gross-Pitaevskii (GP) equation, we refer to 
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the lectures by Burnett in this volume. The condensate density is given by 
n c (r, t) — |<I>(r,i)| 2 and hence (||) reduces to a NLSE for $(r,i), namely 



it) 



at 



2m 



^ + V ex (v)+g\$(r,t)\< 



*(r,t). 



(5) 



Since all atoms are in the identical quantum state, there is no exchange mean 
field in (@). 

The T = GP equation in (JsJ) has been the subject of literally hundreds 
of papers since the discovery of BEC in laser-cooled trapped atomic gases - 
and the equation .appears in almost every chapter in this book. As discussed 
in recent reviews Era, it gives an excellent quantitative description of both the 
static and dynamic (linear and non-linear) behaviour in trapped Bose gases 
below about T < 0.5Tbec- The accuracy can be of the order of a few percent, 
which is as much as one can expect since the non^condensate fraction of atoms 
at T = is also estimated to be about 1% or so. E3 The T — GP equation has 
been extended to deal with two-component Bose gases (involving two different 
atomic hyperfme states) and the effect of perturbations related to laser and rf 
fields (used to manipulate the Bose condensates) is easily incorporated. For 
further discussion and references, see the lectures by Ballagh in this book. 

Our main purpose in this section is to use the GP equation (|j) to discuss 
collective oscillations of a pure condensate. In Section |L we will discuss the 
extension of this equation at finite T to deal with the effect of non-condensate 
atoms. We first briefly review the static equilibrium solution $o(r) of the GP 
equation, given by 



(G) 



where fi is the chemical potential of the condensate. The physics behind this 
can be seen from 

(n - l$(r, t)\N^ = e iEN -^ n (N - 1\^>{t)\n\ e" lENt/n 

= (N - l\y/N\N - Ij e -i(EN-EN-i)t/H 

= y/Ne'^ t/R . (7) 

Using (H) in (||) gives 

h 2 v 2 



^ c0 $o(r) = 



2 m 



+ V ex (r)+g\^ (r)Y 



*o(r). 



(8) 



Assuming a vortex-free ground state, this GP equation for the static conden- 
sate wavefunction < &o( r ) = \J n c o(r) leads to the following expression for the 
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equilibrium condensate chemical potential 

Mco = - - ^-^V + V ex (r) + <?n c0 (r). (9) 
2m^n ca (r) 

A standard approximation in solving (^|) is to ignore the kinetic energy 
associated with the condensate amplitude y/rQ, ie, neglect the — n ^ a term. In 
this "Thomas- Fermi" (TF) approximation, the static GP equation for <I>o( r ) 
reduces to El 

[V ex (r) + g\%(r)\ 2 ] = ,i cQl (10) 
which is easily inverted to give the condensate density profile 

n c o(r) = -[fj, c0 - V ex (r)] 
9 

1 .2„2 



(11) 



_ 1 
9 

Clearly in the TF approximation, the size of the condensate is Rtf, where 

Mco = -mui 2 R^ F . (12) 
One finds /j, c q from the condition j drn c o(r) = N = N c , which gives 

2/5 

; a HO = (h/mto ) 1/2 . (13) 



McO = — 



N c a 
15—^ 



ano 



We recall that the oscillator length ano is the size of the ground state Gaus- 
sian wavefunction of an atom in a parabolic potential. Combining (|l2| ) and 
(EU) gives 



Rtf = a H o 15 



N c a x 



ClHO 

» a HO , if — » 1. (14) 
ano 

The TF approximation ([ll]) for n c o( r ) is very good for large N c , except for a 
small region near the edge of condensate (r ~ Rtf)- In practice, for typical 
values of ano and a, one finds the TF approximation is very good for N c > 10 4 
atoms, a 

To discuss condensate fluctuations around the static equilibrium value of 
$o(r), we first reformulate the time-dependent GP equation in terms of the 
condensate density and phase variables 

*(r,i) = v^M)e i8(r,t) - (15) 



Canberra-rev: submitted to World Scientific on February 1, 2008 5 



Inserting this into the GP equation (||), gives 

ih^-^ - V^c^r = Me(r, t)y/n c - ih 2 ^V 2 + 



dt v L dt ^ ' Jy c 2m 2m K ' 

(Vv^) • V0, (16) 



where we have defined 



h 2 V 



Mc(r,t) = -— *=r- + V ex (r) +gn c (r,t). (17) 

ATTi\i fir 



Separating out the real and imaginary path gives two equations. The real 
part gives 

00 1 
h— = -([i c + -mv 2 ), (18) 

where the condensate velocity field is defined by 

mv c (r,t) = KV9(r,t). (19) 
The imaginary part of ( |l6| ) gives 

i-|^ = -^V-v c -v c VV^, (20) 

which can be rewritten in the form 

= -n c V ■ v c - v c • Vn c = - V ■ (n c v c ). (21) 

^From now on, we shall set Ti = 1 except in some final formulas. 

In summary, we have shown that the time-dependent GP equation for the 
two-component order parameter $(r, t) is completely equivalent to the follow- 
ing two coupled equations for the condensate density n c (r, t) and velocity field 
v c (r,t): 

drlc T7 I \ 

^ + 5^)— <22 » 

where the position and time-dependent generalized condensate chemical po- 
tential n c is defined in (|T7j). In this formulation, a complete description of the 
condensate dynamics is given in terms of two variables n c (r,i) and v c (r,i), 



reminiscent of the hydrodynamic ec 
BEC literature, the equations in 



uations for a classical fluid. In the recent 
are often referred to as the "hydro- 



dynamic" theory. A better description would be to call it the "quantum 
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hydrodynamic" theory, since it is equivalent to the mean-field GP equation. 
In later sections, we shall see that an exact GP equation taking into account 
the dynamics of the non-condensate leads to a generalized set of equations 
quite analogous to (|22|). The key equation ( |l9| ) defining the superfluid ve- 
locity field of the condensate will turn out to be quite general. Needless to 
say, all aspects related to superfluidity of a trapped Bose gas (as compared to 
Bose condensation) are tied to the fact that the condensate exhibits motion 
related to the gradient of a phase (which means that, ignoring vortices, the 
condensate motion is irrotational, V x v c (r, t) = 0). 

Stringari □ first pointed out that within the dynamic TF approximation, 
one could combine the two equations in (|2^) into a single condensate wave 
equation. This approach will be the basis of our development in these lectures. 
Neglecting the kinetic energy term —V 2 ^fh~Z in ([[7]) as small compared to the 
condensate interaction energy gn c , we linearize the resulting equations around 
the equilibrium values 

n c = n c0 + 5n c 

v c = v c0 + 6v c . (23) 



We obtain 



d5n c 

dSv c 
dt 



-V • (n c0 6v c ) - V • (v cQ 6n c ) 

- V (n' cX) + gSn c + mvco ■ Sv e ) , (24) 



where [i' c0 = fj, C Q + |mv^ is independent of position. Assuming that there is 
no vortex in the solution of the static GP equation (ie, v c o = 0), ( pi[ ) reduces 
to the following coupled linearized equations for Sn c and (5v c : 



d5n c 



dt 

38w c g 



V • (n cQ (r)5v c ) 

V5n c . (25) 



in 



ThesfipCan be combined to give the well known T — Stringari wave equa- 
tion m 

^ = iv.Ko(r)Vfo c ]- (26) 

Since the derivation used the TF approximation (where n c o(r) is given 
by ( [ll"|) for r < Rtf), one can rewrite 
isotropic trap potential) 



2q ) in the equivalent form (for an 
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0.5 1 1.5 

N c a/a H0 

Figure 1. Plot of the lowest condensate mode frequencies (monopole and quadrupole) as 
a function of N c . The full curves are the solutions of the coujpled Bogoliubov equations 14 
and the dashed lines arerfehe Stringari TF approximation Ref. B. The data points are for a 
JILA trap with A = v^8liil- For these trap parameters, Nc varies from approximately 10 3 
to 10 4 atoms. For further details, see Fig. 14 of Ref. Q .. 



d 2 Sn c 

dt 2 



[1- 



p2 

JXrp p 



VSn c 



r < Rtf, 



(27) 



where /i c o is given by 
dn c (r,t) = 8n u (r)e- iut o! 



It turns out that the normal mode solutions 
(127]) have frequencies which are independent of the 
interaction strength g or the value of N c . This is a feature of the underlying 

of 



TF approximation, which typically starts to breakdown (as noted 
when N c ^ 10 4 atoms. This is shown by explicit numerical solutions 



earjii 

IB 
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the coupled Bogoliubov equations of motion which describe the normal mode 
solutions of the linearized GP equation when we take the kinetic energy of 
the condensate amplitude ^fn~ c fully into account. For values of iV c <10 4 , 
the normal mode frequencies depend significantly on the magnitude of N c , as 
shown in Fig. |l|. 

We also note that ( p6| ) can be equally well rewritten in terms of the super- 
fluid velocity 5v c or, equivalently, the phase fluctuations.^ This emphasizes 
that the measured condensate density fluctuations E31L3 are directly related 
to the existence of phase fluctuations. The existence of collective modes of a 
pure condensate may thus be viewed already as "evidence" of superfluidity, 
the latter being always a consequence of the phase coherence of the macro- 
scopic wavefunction given in ( |l5| ) which gives rise to the irrotational velocity 



in 



(|19D 



We conclude this section with several examples LU'LrLj of condensate normal 
modes based on solving the T = Stringari equation (|2^) . A wonderful aspect 
about the collective oscillations of a condensate in a trapped gas is you can 
"see" them. As Ketterle has remarked, these condensates are robust - one 
can kick them,, shakp them and these "wisps" of Bose-condensed matter keep 
their integrity. EiE3 The uniform Bose-condensed gas is especially simple, since 
6n c = <5n few e 4 ( k ' r ~" Jt ). This gives 

- A ku) = i^(-k 2 )S nku , (28) 



or uj 2 = c^k 2 , with cr, = y/gn c o/m. This recovers the well-known Bogoliubov 
phonon oscillations □ of a uniform Bose condensate. The neglect of the kinetic 
energy in our TF approximation precludes us frqm obtaining the particle-like 
behaviour at large values of the wavevector k. E3 

The Kohn (or sloshing) mode corresponds to the oscillation of the centre- 
of-mass of the static condensate profile with the trap frequency ujq. This mode 
is described by 



n c {r,t) = n c0 {r-r,(t)), (29) 



where = v c (t) and 



d 2 V (t) 



w 2 V {t). (30) 



dt 2 

The proof is simple. Linearizing (p9|), we find 

ra c (r, t) = n c0 (r) - r] ■ Vn c0 (r), (31) 
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which gives the explicit form for the condensate fluctuation, 

SnJr, t) = —muskx ■ rf(t). 
9 

Inserting this into the Stringari equation (| 

d 2 r,(i) 



(32) 



Vn c0 • 



we obtain 



i c0 V( V (t) ■ r)] 



(33) 



dt 2 ^ 

= ujI V • [n cQ r]} 

= ulri ■ Vn c0 . 

This confirms that the centre-of-mass position r)(t) of the static condensate 
distribution satisfies the SHO equation ( |30| ) with frequency ojq. 

The breathing (or monopole) condensate normal mode corresponds to a 
velocity fluctuation of the form 

6v c (r, t) = Ay e~ iut (r < R TF ). (34) 

Using V • r = 3 and V • (r 2 r) = 5r 2 , it is easy to verify from the continuity 
equation in (|2^) that the associated density fluctuation is 

- iuj5n c = -V • [n c0 (r)(5v c ] 



9 



A 



3-5 



B 2 

T F 



or 



<5n w (r) = B 



3 Hjrp 



(35) 



(36) 



where /i c o is given by (|l2]). Inserting this into the Stringari wave equation 



nvcs 



w 2 <5n w (r) 



p2 



5 2r 



3 f^'-p jp 



-5w 2 (5n w (r). 



(37) 



Thus the breathing mode for an isotropic parabolic trap has a frequency lo = 
^/Euiq. As noted earlier, this mode frequency is not explicitly dependent on 
the interaction strength g. However the underlying theory is very dependent 
on mean-field effects. Moreover, we recall that a non-interacting trapped Rose 
gas has a breathing mode with a frequency u> = 2u>q at all temperatures. □ 

As a final example, we consider the so-called "surface" modes of a T = 
condensate, as described by U 



5v c (r,t)=AV [r l Y lr 



(38) 
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This mode corresponds to phase fluctuations of the form 69^ = 
mAr l Yi m (9,4>)- One may easily verify that V • i5v c = and hence from 
(Eq) we obtain 



This is equivalent to 



d 2 5v c 
dt 2 

d 2 69 



wgV(Jv c -r). (39) 



-u4{V86)-y (40) 



or 



-oj 2 5d u (r) = -u%160 u (t). (41) 

Thus the surface oscillations of the condensate phase have a frequency given 
by uj — VIluq(1 = 1, 2, 3, . . .). 

One great advantage of the quantum hydrodynamic formalism (within 
the TF approximation) is that it is easy to also treat the T = normal 
modes of an..apisotropic parabolic well. We refer to the literature for further 
discussion. BffiLj We also note that it is straightforward to derive a Stringari 
wave equation for the fascinations of a vortex state (this has been done by 
Svidzinsky and Fetter Ej). 



3 Generalized GP equation at finite temperatures 

In this section, we generalize the T — GP equation of Section || to finite 
T, so that it includes the effect of the coupling of the condensate to the 
non-condensate degrees of freedom. We concentrate in this Section on how 
the T = condensate modes are renormalizcd and damped by the thermal 
cloud. We defer discussion of the appearance of new collective modes mainly 
associated with the thermal cloud to sections ^ and ||. 

Up to the beginning of 2000, the study of the dynamics of a trapped Bosc 
gas at finite temperatures has been largely ignored by experimentalists but 
actively studied by many theorists (especially in the last year or so). One 
of the reasons for the lack of finite temperature data is that there are so 
many interesting phenomena to study at T ~ 0! However, another reason 
seems to be the implicit belief that the presence of the non-condensate just 
complicates the behaviour of a pure T = condensate - but is not the source 
of any interesting new physics. In this and subsequent sections, I hope to 
argue that this "belief" is quite wrong. The coupling of the condensate and 
non-condensate degrees of freedom at finite T leads to a new two-component 
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system in which both components can exhibit coherent behaviour, analogous 
to the well-known superfluid macroscopic phenomena appearing in liquid 4 He. 

We shall see that, as expected, the finite temperature GP equation of mo- 
tion for Q(r,t) is not closed. It's general solution involves knowing the equa- 
tions of motion for the non-condensate. We will work within an approximation 
where the non-condensate atoms can be described by a quantum Boltzmann 
transport equation for the single-particle distribution function /(p, r, t), with 



"( r >*) = / 7^W/(P>M)- ( 42 ) 



dp 

(2n) 3 

The details of this kinetic equation will be developed in Section [| In this 
section, we concentrate on incorporating the effect of a static thermal cloud 
into the finite T dynamics of the condensate. 

The theory of interacting Bose-condensed fluids is most usefully discussed 
using quantum field operators. This .procedure was formalized-hy Beliaevtj in 
1957 and developed by Bogoliubov E_Zl, Gavoret and Nozieres Ea, Martin and 
Hohenberg Ej, and others in the early 1960's. We recall: 

tp' (r) = creates atom at r 

■0(r) = destroys atom at r. (43) 

These quantum field operators satisfy the usual Bose commutation relations, 
such as 

ip(r),ft(r')\ =S(r-r'). (44) 

All observables can be written in terms of these quantum field operators, such 
as the density n(r) — ip' 1 (r)ip(r) and interaction energy 

= ~<? / dr^( r )^t (r) ^ (r)v 3 (r) . (45) 



The crucial idea due to Bogoliubov caO and later generalized by Beliaev U 
is to separate out the condensate component of the field operators, 

$(r) = (^(r))+^(r), (46) 

where 

■0(r)) = $(r) = Bose macroscopic wavefunction. 
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This quantity plays the role of the "order parameter" for the Bose superfluid 
phase transition: 



$(r) =0 T>T C 
^ T <T C . 

We note that $(r) = *Jn~ c e lB is a 2-component order parameter. Clearly, 
3>(r) is not simply related to the many-particle wavefunctions \t(ri,r2, . . . rjy). 
The thermal average in /^(r)^ involves a small symmetry-breaking pertur- 
bation to allow <I>(r) to be finite, 



Hsb = lim / dr 

ri— >0 



r)(r)ft(r) + r)*(r)4>(r) . (47) 



The philosophy behind the concept of symmetry-breaking was extensively 
discussed by Bogoliubov EJ, for a variety of condensed matter systems, in an 
article which is still highly recommended. 

It is useful to make a few comments on the physics behind $(r, t). $(r, t) 
is a coherent state, with a "clamped" value of phase - rather than a Fock- 
state of fixed N, with no well-defined phase. <£>(r, t) acts like a classical field, 
since quantum fluctuations are negligible when N c is large. Probably P.W. 
Anderson Ej deserves the greatest credit for understanding (in the period 1958- 
1963) the new physics behind working with a broken-symmetry state $(r,i), 
both in BCS superconductors and in superfluid He. It captures the physics 
of the new phase of matter (such as the occurrence of the Josephson effect) 
and the associated superfluidity. The symmetry-breaking perturbation ( [47j ) 
allows the system to internally set up off-diagonal symmetry-breaking fields, 
which persist even when the external symmetry-breaking perturbation in ( ff7j ) 
is set to zero at the end (?7 — ► 0). The same sort of physics is the basis of the 
BCS theory of superconductors. 

One can formulate the GP and Bogoliubov approximations directly in 
terms of a variational many-particle wavefunction (see the lectures by Leggett 
in this volume). However, such formulations are limited to simple mean- field 
approximations. The explicit introdiic±km_o£_thc broken-symmetry order pa- 
rameter $(r, t) gives a systematic way BOOEi of isolating the role of the Bose 
condensate in a general treatment of an interacting Bose-condensed fluid. As 
we shall see, the resulting formalism allows one to deal with questions related 
to damping as well as superfluidity in both the collisionless and hydrodynamic 
regions. 
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The exact Heisenberg equation of motion for the field operator is 



dip(r, t) 

l ^t~ 



2m 



+ V ex (r)+SU(r, t) 



V>M) 



+ 7 ? (r)+ 3? /- t (r,t)V3(r,i)^(r,i), (48) 
where SU (r, t) is a small time-dependent driving potential. This gives an exact 
equation of motion for <&(r, t) = ($(r, t) \ , 



.<9$(r, t) 



2m 



+ V ex (r) + SU(v, t) 



*(r,t) 



+ V( T ) + 9 (^ f (r, t)-0(r, t)-0(r, t) 

with [(using the decomposition (phf)] 

$t$^ = |$|2$ + 2 |$| 2 ^ + $ 2 ^,t + + 2$^t^ + ^t^_ 

Taking the symmetry-breaking average of ([30]), one finds 
where 

n c (r,t) = |$(r, i)| 2 = condensate density 
n(r, i) = (ffi(r, t)-0(r, i)^ = non-condensate density 

m(r, i) = (i[>(r, i)V>(r, t)^ = off-diagonal (anomalous) density. 

Using (Ell) in (fi9|), our "exact" equation of motion for <I>(r,i) is 



(49) 
(50) 
(51) 



(52) 



<9$ 



2m 



+ V ex + gn c (r, t) + 2gh(r, t) 



+ gm(r, t)<P* + g (ft(r, t)^(r, t)i>(r, t) 



(53) 

We now consider various approximations to the generalized GP equation in 
(1): 

(a) The Hartree-Fock-Bogoliubov (HFB) approximation for $ corresponds 

to neglecting the three-field correlation function (ift but keeping 

the n c , h and m fluctuations. This HFB has been exhaustively treated in 
many recent papers- in conjunction with the separate equations of motion 
for n and m. This HFB approximation can be used to generate the 

same normal mode spectrum which the Beliaev second-order self-energy 
approximation gives at finite T. cB 
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(b) The dynamic Popov approximation corresponds to ignoring both (^'ip^jjj 

and m — (jP^P) m (HI)- Theories of this kind involve coupled equations 
for $ and n. c3 



(c) The static Popov approximation involves even a further simplifica- 

tion, namely it ignores fluctuations in the density n(r,t) of the thermal 
cloud, 



n(r,t) ~ n c0 (r). 



(54) 



Using in (^3|) corresponds to treating the dynamics of the condensate 
moving in the static mean field of the non-condensate thermal cloud, ie, 



ih 



<9$ 
~dt 



2m 



+ V ex (r) + 2gh (r) + gn c (r, t) 



®(r,t). 



(55) 



There is a considerable literature based on this static Popov-approximation. 
In a nutshell, in this section, we will discuss an extension of fl55| ) which involves 
including the damping associated with the {y 1 term, but again treating 
the thermal cloud statically. 

As in Section |^, we first rewrite ( |53| ) in terms of the condensate amplitude 
and phase. Inserting ( [L5| ) into d5^) and following the same procedure which 
led to the T = equations in (|22|), we obtain: 



~dt 



V • n c v c = 2glm 



$* 2 m + $*( t /i t #) 



(56) 



80 
di 



(V6 



2m 



(57) 



where now (compare with (jl^)) 



Mc(r,t) 



2m^fn 



+ Vex + gn c + 2gn 
$* 2 m + fffi<ijnp\ 



(58) 



Formally, (|53| ) and hence ([56|) - (58) are exoci results. In these lectures, we will 
limit ourselves to finite temperatures where the dominant thermal excitations 
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are well approximated by high energy non-condensate atoms moving in a self- 
consistent dynamic HF mean-field 

P 2 

i p (r,t) = £ r +I4*(r)+2#[n c (r,i)+n(r,t)] 



2m 

2 

P 



2m 



(59) 



In this region, we shall neglect the anomalous pair correlations among the 
thermal atoms (ie, fh = 0). One can show that both m and (^ipip) in (56) 
and (|5^) are of order g. Thus one sees the correction terms in ( p6| ) and (58) 
involving these functions are of order g 2 . In fact, calculation shows that to 
order g, (ip^ ipip) is imaginary. Our procedure is to ignore the 0(g 2 ) terms 
in ([36]) and (|5^) except those which are imaginary and hence are a source of 
damping of condensate motion. 

The end result of this approach is that we are left with 

V-n c v c = -r 12 [/,$] 



with 



Mc(r,*) 



Of 
dv c 
dt 

^ 2 \/n~c 



2m^/n, 

and we have defined the new function 

T 12 [f,^} = -2glm 



+ V ex (r) + gn c (r, t) + 2gh{r, t) 



(60) 



(61) 



(62) 



We see that Ti2 plays the role of a source term in the condensate con- 
tinuity equation in (|60| ) . It depends on both <5(r, t) and the single-particle 
distribution function /(p,r, t). It is useful to note that and ( |6l|) are 
equivalent to the approximate GP equation (see also Refs. oEj) 



if i 



h 2 v 2 

2m 



+ V ex {r) + gn c (r, t) + 2gn(r, t) - ihR(r, t) 



where the dissipative function is 
R{v,t) = 



ri 2 [/,$] 

2n c {v,t) 



0(g 2 



(63) 



(64) 



To proceed, we need to calculate the field correlation (-0t-0V>) n which 
termines F12 in (|6^). This has been done in Appendix A of Ref. 13, working 
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to lowest order in g at finite T where the single-particle spectrum in (59) is 
adequate. The result is 



(^*4nP) = -ig ^0L J d Pl J dp 2 J dp 3 



x S (p c + pi - p 2 - P3)<5(e c + h - e 2 - e 3 ) 

x [/i(l + / 2 )(1 + h) - (1 + h)hh] , (65) 

where /i = /(p, r,i),e c = fi c (r,t) + |mv^(r, i) is the local condensate atom 
energy and p c = rnv c is the condensate atom momentum. In Section EL we 
shall see that ([35]) is closely related to the collision integral C\ 2 [/, $] which 
enters the Boltzmann equation for /(p, r, t) and describes collisions with the 
non-condensate atoms in which one condensate atom is involved. We note 
that to leading order, the expression in ( |65| ) is imaginary. Using this in (p^), 
we find 

r 12 (r, t) = 2ff 2! ^^ J dpi J dp 2 J dp 3 S(p c + pi - p 2 - p 3 ) 

x 6 (e c + ei - e 2 - e 3 ) 

x [A(l + / a )(l + / a ) - (1 + /O/2/3] • (66) 

As a first application of our new equations in (|60|), we will ignore the 
dynamics of the non-condensate cloud. That is, we will assume that the 
condensate interacts with a static thermal cloud in thermal equilibrium 



/(p,r,t) ~ /°(p,r) = e0[p2/2m+UoM _ M _ 1 , (67) 

where p,Q is the equilibrium chemical potential of the non-condensate atoms. 
In Section ^, we will discuss the C22 collisions which produce this static equi- 
librium Bose distribution. The value of /to is known since it must equal the 
condensate equilibrium chemical potential fj, C Q given by the static solution of 
the generalized GP equation. Within the TF approximation, mWj reduces to 



Mco = V cx (r) + gn c0 (r) + 2gh (r) (68) 
and hence we see that (|59| ) reduces to 

2 

P 

e P ( r ) - Ao = h U (r) - Jxq 

= ^+gn c0 (r). (69) 
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This result is consistent with the correct Bogoliubov "excitation energy." lilli 
The static non-condensate density in ( |68| ) is given by 

no(r) = ^.g 3/2 (*o = e -/Wr)) . ( 70 ) 

We can now calculate ri2(r,i) using the equilibrium Bose distribution 
( |67| ) for the thermal cloud. One finds B (for a related calculation, see Ref. El). 

r° 2 (r,i) = ^||[ e -fl(fio-«.('.*)-i™J) - i], (71) 

where we have defined a C\i collision time 
1 _ 2 3 2 



Jdp! Jdp 2 Jdp 3 S(p c + pi p 2 p 3 ) 



ri 2 (r,t) (27r) 5 , 

x 5(e c + e pi - e P2 - e P3 ) (1 + /fj/J/a ■ (72) 



ri2[/o, < J>] in ( |7l| ) still depends on n c (r, i) and v c (r,i) of the condensate 
through fi c , p c and e c . We note that when the condensate is in static equilib- 
rium, /i c (r, t) — * /i c o and v c o = 0. In this case, /i c o = fla and the expression 
in the square bracket in ( |7l| ) is seen to vanish. Thus r 12 (/o, ( i > o) = 0, as it 
should when both components are in static thermal equilibrium. 

Summarizing, at this point we have a closed set of equations which can 
be used to describe the dynamics of the condensate in a trapped Bose gas 
which include the interactions with a static equilibrium thermal cloud. These 
equations are 

^ + V.n c v c = -r 12 [/ ,$] (73) 

m (J^ + v c • Vj v c = -V/x c , (74) 
where T^/ , $] is given explicitly by the expression in ( |7l| ) and 

Hc(r,t) = - 2m V + V ex (r) + 2gn (r) + gn c (r,t). (75) 



^From ( |73[ ) and (|74|), we can obtain linearized equations of motion for the 
condensate fluctuations 8n c and <5v c . We use the fact that, to lowest order in 
the fluctuations from static equilibrium, ([FT]) reduces to 

dT 2 = M£ 6flc{r>th (76) 

T 12\ r ) 
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where the "equilibrium" C\i collision rate [using (|69|) in (|72[)] is defined by 
1 _ 2 ff 2 



,(r) " (27T)' 1 



J dpi J dp 2 J dp 3 5(pi - p 2 - pa) 



x 5{ Pl 2 1 P§ ~ ^co)(l + /i )/ 2 /3 • (77) 

In the static TF approximation, we recall that the equilibrium distribution is 
given by /? = [e^O* / 2m +s™oo) _ 

In the present discussion, we further restrict ourselves to the dynamic 
Thomas- Fermi limit valid for large N c , 

+ V • (n c0 5v c ) = -^Sn c (78) 

TO ^T = ~ gV6ric ■ ^ 

The collision time r'(r) describes collisions between the condensate and non- 
condensate atoms when the condensate is perturbed away from equilibrium, 

1 _ g"rf>( r ) 1 rsn x 

r'(r) fc B T 7f a (r)- l8Uj 

The new term on the right-hand side of (|7^) causes damping of the conden- 
sate fluctuations due to the lack of collisional detailed balance between the 
condensate and the static thermal cloud. We note that this collision time 
t 1 (r) is only a function of the position r through its dependence on the static 
condensate density n c o(r). 

We can easily combine (f78|) andJ79|) to obtain what we shall refer to as 
the finite T Stringari wave equation a 

d 2 5n c g 1 dSn c 

-— V • (n c0 V<5rc c ) = — 7 — — . 81 

ot z m t at 

If we neglect the right-hand side, we obtain the undamped finite T Stringari 
normal modes Sn c (r,t) = Srii(r)e~' luJit given by the solution of 

- -V • [n c0 (r)V^(r)] = ^Sn t (r). (82) 

TO L J 

As has been noted by several authors in recent papers EH, n c o(r) at 
finite T can be well approximated by the TF condensate profile at T = 
but with the number of atoms in the condensate N C (T) now being a function 
of temperature. This is because the static mean field of the non-condensate 
plays a minor role. In the regions where the condensate density is finite, we 
effectively always have n c o(r) 3> n-o(r), as illustrated in Fig. @. With this 
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2 4 6 8 
r(a HO ) 

Figure 2. Density of atoms as a function of position in an isotropic parabolic trap containing 
5000 87 Rb atoms at T = 100 nK, with T BE c = 14 9 nK. At this temperature, the total 
number of atoms in the condensate and non-condensate are equal. The solid line is the 
total density and the dashed line is the non-condensate density j-i The TF approximation is 
not used. The parameter ano is the oscillator length. See Ref. 0. 



approximation for n c o (r) , the solutions of the finite T Stringari equation ( pl| ) 
will be identical to those at T = (see Section |^) , since the T = Stringari 
frequencies do not depend on the magnitude of N c . OLjCOurse, as shown 
by calculations solving the coupled Bogoliubov equations EniJ, when N c < 10 4 
the TF approximation breaks down. Thus the condensate collective mode 
frequencies will always become temperature dependent close to Tbec> where 
the TF approximation is no longer valid. E3 
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We can use the undamped Stringari modes given by ( |S2| ) as a basis set to 
solve ( |8l| ) and find the damping of these modes. Writing <5n c (r) = ^ c,£rii(r), 
and using the orthonormality condition f drSni(r)Srij(r) — 5ij, one obtains 
the following algebraic equations for the coefficients 

u?c t = w?Cj - 7yCj , (83) 

3 

where 

7 y= / ' dT8ni{v)6nj{T)/¥{T). (84) 
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Assuming the damping is small (our present discussion is for the col- 
lisionless region), ( |S3| ) is easily solved using perturbation theory by setting 
7y = for j 7^ i. This gives the damped Stringari frequency (to lowest order) 
f2j = u>i — iTi, with 

This result for i\ is reasonable, namely it involves a spatial average over 
l/r'(r) weighted with respect to the undamped density fluctuations of the 



Stringari wave equation (81). Calculation shows that the effect of coupling to 
other modes ("yy ^= 0) is extremely small. 

Before discussing the trapped gas, it is useful to first apply our theory 
to a homogeneous gas, where r' is independent of position. In this case, 
( |S5| ) simply reduces to Tj = l/2r'. Although our model applies only to the 
collisionless region, it is useful to compare the inter-component collision time 
for a uniform gas in the collisionless and hydrodynamic regimes. In Section 
we show that the inter-component collision time in the hydrodynamic 
region El is given by r M = ujjt', where the temperature-dependent factor 
<jh depends on various thermodynamic functions [cry is given explicitly by 



(165) in Section [4j. In Fig. |3|, we compare \jr' and 1/t m as functions of T. 
We see that cx# dramatically alters the inter-component relaxation rate 1 /r„ 
appropriate to the hydrodynamic regime, as compared to 1/r' involved in the 
collisionless regime. For completeness, in Fig. p| w e also plot the often-used 
classical collision time as well as rf 2 defined in (|77|) 

We briefly discuss some numerical calculations □ of the inter-component 
damping using the above formalism. We consider 87 Rb atoms in a spher- 
ically symmetric trap with frequency vq = 10 Hz and N = 2 x 10 6 . In 
the collisionless limit, we require 0JiT c \ 3> 1. We obtain an upper limit on 
1/t c i by taking the density in the center of the trap n(r — 0), which gives 
1/t c i = &a 2 N u^m / {irksT) . For the parameters we use, lo\qt c \ 19 (coat 
pared to uj m T c \ ~ 20 for the JILA dataEJ and oj 2T c \ ~ 2 for the MIT data El 
on collective oscillations at finite T). 

In Fig. U we plot the damping rate r 10 for the breathing mode (n = 1, 
I = 0) as a function of temperature up to T = 0.95TbeCi where N c ps 7 x 10 4 . 
At higher temperatures, the Thomas-Fermi approximation will start to break 
down and the mode frequencies become temperature dependent. 

The damping of condensate modes we have discussed in this section is 
due to the fact that the condensate is out of diffusive equilibrium with the 
thermal cloud. At finite T, the collective oscillations of the coupled conden- 
sate and non-condensate can be generally split into two classes. One mode 
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Figure 4. Breathing mode damping rates vs. temperature. In these plots, the damping rates 
are normalized to the mode frequencies and we only plot up to T = 0.95Tbec> above which 
the Thomas-Fermi approximation will start to break down. We show the inter-component 
damping rate given by ( |8o| ) and (pj) for the breathing mode (n = 1, 1 = 0) of a spherically 
symmetric trap, as well as the Landau damping. Hj For further details, see Ref. Q. 



mainly involves (out-of-phase) motion of the condensate, and is the finite T 
generalization of the analogous T = condensate mode, cl For the same sym- 
metry, there is another mode which mainly involves the (in-phase) motion 
of the non-condensate thermal cloud. This mode is naturally viewed as a 
generalization of the T > Tbec thermal cloud oscillation to temperatures 
below jCbbc- These two types of modes (for a given symmetry) have been ob- 
tained 13 in the collision-dominated hydrodynamic region, when the dynamics 
of the thermal cloud is fully allowed for. Treating (to first approximation) 
the thermal cloud (the non-condensate component) as always in static equi- 
librium in this section means that our theory is only applicable to condensate 
oscillations which involve motions which are out-of-phase with the thermal 
cloud. 

As we discuss in Section ||, the possibility that the condensate may be out 
of diffusive equilibrium with the non-condensate is a general feature of the 



dynamics of trapped Bose-condensed gases, 
region, it leads to a new relaxational mode. 



the hydrodynamic two-fluid 
This effect is missed in the 
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well-known form of the two-fluid hydrodynamics developed by Landau &Lj, 
where one (implicitly) assumes that the normal and superfluid components 
are always in local diffusive equilibrium. We will return to this question in 
Sections [I] and ||[ 

Of course, in addition to the inter-component relaxation discussed in 
this sectian-i-flpe also has Landau and Beliaev damping of condensate os- 
cillations. ESEjGj These mechanisms are treated in more detail by Burnett in 
this book. Landau damping is briefly discussed in Section ^. 

Clearly one can extend the formalism of this section to any problem based 
on the T = GP equation. In particular, it can be used to discuss this kind 
of inter-component damping at finite T of oscillations in two-component Bose 
gases, vortex dynamics EjEj, and Josephson oscillations between two traps. 



4 Dynamics of the coupled condensate and non-condensate at 
finite temperatures 

In contrast to Section we now grapple with the dynamics of the thermal 
cloud. As we noted in Section in dealing with this very complicated prob- 
lem, we will be quite modest and treat the non-condensate using the simplest 
microscopic model which captures the important physics. We limit ourselves 
to finite T, where the non-condensate atoms can be described by the particle- 
like HF spectrum given by (p9|). For trapped Bose gases, this spectrum is 
probably adequate down to quite low temperatures, for reasons discussed in 
an important paper by the Trento group. Lj To extend our present analysis to 
very low temperatures is much more involved since then the excitations of the 
thermal cloud take on a collective aspect (ie, a Bogoliubov-type quasiparticle 
spectrum must be used). Ej 

The single-particle spectrum ( |59| ) is appropriate in the semi-classical limit, 
where we can use the single-particle distribution function /(p,r,i) given by 
the solution of a kinetic equation. This procedure generalizes the approach 
of Boltzmann (1880's) for describing binary collisions in a classical gas. Such 
a quantum Boltzmann equation for a trapped Bose-condensed gas at finite 
temperatures has. been derived and extensively discussed by Zaremba, Nikuni 
and the author. 13 The conditions of validity are 

ksT > gn c , fc^T > hui , (86) 

where Tiloq is the spacing of the energy levels of the harmonic trap. M<ffl3 
general but less explicit discussions of kinetic-equations are given in Refs. EJcj. 
Related work of Gardiner and coworkers cB is discussed in the lectures by 
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Ballagh in this book. The quantum kinetic equation we use is given by 

WPllA + £ . V r /(p,r,i) - V r U(r,t) ■ V p /(p,r,i) 
at m 

= C 22 [f]+C 12 [f}. (87) 

The right hand side describes how binary collisions effect the value of the 
single-particle distribution function /(p, r, t). The effective time-dependent 
HF potential U (r, t) is defined in (^) . The effect of collisions between excited 
atoms in the non-condensate is described by: 

C 22 [f] = J dp 2 J dp 3 J dp 4 5(P + p 2 p 3 p 4 ) 

X 5 (ip + £p 2 — £p 3 — £p 4 ) 

x [(1 + /)(1 + f 2 )hh - ff 2 (l + / 3 )(1 + U)] . (88) 

We recall that creating a Boson gives a factor (1 + /) and destroying a Boson 
gives /. In the classical high temperature limit, / < 1 and the collision integral 
C 22 in ( |8"s| ) becomes much simpler. 

In addition to C 22 collisions, we also have collisions which involve one 
condensate atom: 

Ci 2 [/] = yf^ja J dpi J dp2 J dp^i^c + Vi -P2 -P3) 

x S (e c + £ p i - £p 2 - i p3 ) [5 (p - pi) - S(p - p 2 ) - <5(p - p 3 )] 
x K(l + h)hh - n c /i(l + / 2 )(1 + /s)] • (89) 
For convenience, we recall that the condensate atom has 

1 2 

energy: e c = \i c + -mv c ; fi c = V ex + gn c + 2gn 
momentum: p c = mv c (90) 
We note the key difference between C\ 2 and C 22 collisions: 

• C22 and C\ 2 both conserve energy and momentum in collisions. 

• C\ 2 does not (but C22 does) conserve the number of condensate atoms. 
C12 describes how atoms are "kicked" in and out of the condensate. 



It turns out in the generalized GP equation given in (|63[), the damping 
term —iR(r,t) is closely related to C\ 2 [f]. This makes sense, since the C% 2 
collisions modify the condensate described by $(r, t). One easily may verify 
that [comparing ( |66|) and 



r 12 [/,*]=y Ac 12 [/(p > r > t)]. (91) 
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More precisely, the three-field correlation function in the exact equation of 
motions [see (|5^) and (|58|)] is related to C\2 [/] by 

/ j§r Cl2[f] = 2gVTFc Im (^^) • (92) 

We have to solve the equation of motions for /(p, r, t) and <i>(r, t) treating 
C12 [/] very carefully. There will be an exchange of atoms between the n(r, t) 
and n c (r, t) components through the C\i collisions. We can use these coupled 
equations for a variety of problems. Recently these coupled equations have 
been used in two problems: 

(a) Derivation of a generalized set of two-fluid hydrodynamic equations. @ 

(b) Discussion of the rate of growth of a condensate @ due to a sudden 
quench in which the high energy spectrum of the thermal cloud distribu- 
tion is suddenly removed. This work has many points of contact with-thc 
formalism reviewed by Ballagh in this book for condensate growth. cB 

In this section, we will linearize these equations and consider the collec- 
tive oscillations of the combined system composed of condensate and non- 
condensate. l3 It is usefullo introduce two regimes to describe collective modes 
in interacting systems ElO: 

I. Collisionless (produced by mean fields) 

or T < t r (uj = (93) 

II. Hydrodynamic (produced by collisions) 



uit r < 1 or T > t r , (94) 

where tr is some appropriate relaxation time. What .should we use for tr? 
For a uniform classical gas, this is the collision time E3'E3 

1 = V2nav (95) 



where 



a = 87ra 2 (for Bose atoms); a = s-wave scattering length. 



v ~ average velocity of atoms ~ \ — — . 

V m 

h = density of excited atoms. 
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Even for a Bose-condensed gas, taking tr ~ t c is a reasonable first estimate. 
To get into the interesting hydrodynamic region (ljtr <C 1), we need a small 
value of tr, ie, a large density n or a large collision cross-section a (perhaps by 
working near a Feshbach resonance) . We note that current BEC experiments 
deal with N ~ 10 6 atoms or larger and thus the kind of density of thermal 
atoms needed to be in the hydrodynamic region now seem achievable. 

Let us look at the kinetic equation (p7h, writing it in the schematic form: 



In the collisionless region, we need only solve Cf = 0. For comparison with the 
collision-dominated hydrodynamic domain, we briefly discuss some features 
of this collisionless domain, limiting ourselves to a uniform gas for simplicity. 
This enables us to understand how Landau and Beliaev damping arise from 
the collisionless Boltzmann equation, in contrast to the inter-component col- 
lisional damping (discussed in Section |J) which has its origin in the C\i [/, $] 
collision integral. __ . 

We consider linear response theory EZrE3 for the deviations of /(p, r, t) 
from thermal equilibrium 



For a uniform normal Bose gas to first order in Sf, the collisionless kinetic 
equation reduces to 



—L + P. . y rS f _ V r [2gSn(r, t)} + SU(r, t)} ■ V p / (p) = 0. (98) 
at m 

Assuming SU(r,t) = 8Uk^e l( - k ' r ~ tAjt \ then we have £/(p,r } i) = 
<5/fcw(p)e i(k r "" t) , where Sfku(p) satisfies 




(96) 



/(p,r 1 t) = / (p J r) + 5/(p,r,t). 



(97) 



iu)Sf(p) + i 



pk 



6f(p) - [2gSn kLU + SU kL ,}ik ■ V p / (p) = , 



(99) 



m 



where 



/o(p) 



1 



/o(ep) 



(100) 



and 




(101) 



Solving (gg), one finds 




(102) 
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or 



Here 



Sn k« = X l m ^"l . <5t/ fccJ = Xnn(k,Lu)5U kuJ . (103) 



is the fc — > limit of the density response function of a non-interacting uniform 
Bose gas. 

This shows how one can use the collisionless kinetic equation to find the 
density response function Xnn(k, uj) of a unif orm Bose gas. The density fluctu- 
ations are given by the poles of x„„(k, u>) in (103). In this RPA with exchange, 
the solutions of 

l-2 5X ° n (k,u,)=0 (105) 

are called "zero sound" modes (following Landau's terminology in the analo- 
gous calculation for an interacting Fermi gas EJ) . It is clear that the solutions 
u){k) of (105) will be Landau damped, with a width related to 

- 2< ? Jmx° n (k,w + iO + ) 
dp 
(2n) a 

Clearly this damping will come from the zero sound mode absorbing a thermal 
excitation of energy e p to create a thermal excitation e p +fe, 

huj(k) = e p+ k - e P •• (107) 

We might also remark that if this linear response calculation is extended 
to deal with a Bose-condensed gas, one finds a new damping mechanism be- 
cause now XnnO*-^) describes non-interacting Bogoliubov excitations E p . In 
addition to the contributions of the kind given by (107), one finds x^ n (k, u>) 
now has additional poles at 

hu(k) = E p+k + E p . (108) 

These give rise to "Beliaev" damping, where the zero sound mode can decay 
into two Bogoliubov excitations. In the limit of T — ► 0, Beliaev damping 
dominates over Landau damping. It was first calculated at T — by BeliaevllJ 
in 1957 in his classic study of a weakly interacting Bose gas. 

Several recent papers have discussed. Landau and Beliaev damping of col- 
lective modes in a trapped Bose gas. ElOO'Cj'C3 This problem has been nicely 



2 9 1" / T^l/ofe) - fo(e P +k)]S[ui(k) - (e p+k - ep)]. (106) 
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formulated directly in terms of response functions by Minguzzi and Tosi El in 
conjunction with mean fields produced by the condensate and non-condensate 
density fluctuations. Giorgini l3 has given a physically appealing discus- 
sion working directly with the coupled mean field equations of motion for 
<5(r, t),n(r, t) and m(r, t). The approach of the Oxford group is reviewed in 
Burnett's lectures in this book. 

In this connection, a simple treatment of Landau damping of condensate 
collective modes can be given starting from our generalized GP equation in 
(|63|). In Section pi we ignored the fluctuations in the HF term 2gh(r,t) and 
concentrated on the damping associated with the — iR{r, t) term in (|69|). A 
simple way of including the 5n(r, t) fluctuations induced by tha_mean field 
associated with the condensate fluctuations is to use (see Refs. OH) 

5h(r, u) = J dr'x° m (T, r', to)2g6n c (r', to). (109) 

Here Xfth ^ s ^ ne density response function for a non-interacting gas of atoms 
with the equilibrium HF spectrum given by (|6^) and the chemical potential 
/[to is equal to [i c o in (|6^). For a uniform Bose gas, the condensate collective 
modes given by (B3h are found to be 



to 



2 _ „2 ; 2 



i to 



c ^[i + 4 ffX y fi (k, w )]-- . (no) 



Here cq = ^/gn c o/m is the Bogoliubov phonon velocity Jseii (|28|)1 and r' is 
defined in (pO|). In the long wavelength limit, one finds ffla 

lim -Imxln(k,to = c fc) = (111) 

where the 1/3 factor arises from the exchange term in the single-particle self- 
energy. Thus one obtains 

to = c k - i{T L + -L), (112) 

where the Landau damping is given by (RPA with exchange) 

V L = ^ak B Tk. (113) 

The exact r"^'^ for T l obtained from a careful evaluation of the Beliaev self- 
energies ESejcj at finite T is the same as (|H3|), apart from a slightly smaller 



numerical coefficient (4/3 is replaced by 37r/8). This difference is due to our 
neglect of thp_ffuctuations associated with the anomalous pair density m(r, t). 
Such terms Ej are in the exact equation for $(r,t) given in but were 
neglected in deriving (|63|). 
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To conclude this discussion of different kinds of damping in the collision- 
less limit, we can also include the effect of the collision terms in (87) as a 
perturbative correction to the solution determined by Cf = 0. For a Bose_gas 
above Tbec, one might use the simple "relaxation time" approximation 



The inclusion of such a term in (|9 

dSf 
dt 



C 22 [/]^-(/-/o)/t 22 = -^/t 22 . 

simply involves the change 
Sf 



m 

dt 



T~22 



(114) 



(115) 



and hence one fi nds ( 103 ) again, with the replacement u> — ► u + i/r 22 in the 
denominator of ( 104|). The collision time t 22 is the analogue of the classical 



expression given in ( |95| ) . One sees that including C 22 [/] in the simple relax- 
ation time approximation (114) only slightly alters the Landau damping given 
by ( |106|) . This is quite different from the contribution arising from Ci 2 [/] in 
so far as this enters directly into the generalized GP equation (see Section ||) . 

We now leave the discussion of the collisionless domain and turn to the 
collision-dominated hydrodynamic region. In this case, the collisions between 
the non-condensate atom are assumed to be sufficiently rapid that the form 
of /(p, r, t) is largely determined by the requirement that 



C 22 [/, $] =0. 



(116) 



That is, the collisions are so strong they force the system to be in local equilib- 
rium. The unique solution (denoted by /) of the equation (116) is well-known 
to be given by 



/(P,r,i) 



1 



^- + [/(r,t)-AO,t)] 



1 



(117) 



where v„ is the average local velocity and ft is the local chemical potential of 
the thermal atoms. This local equilibrium Bose distribution involves the local 
variables 0,v n ,ft and U, all of which depend on (r, t). Why must / have the 
form in (117)? To satisfy C 22 [/i] = 0, we must have [see 



(1 + /i)(l + h)hh - + /s)(l + h) = 0, 



(118) 



and this requires that / be given by the Bose distribution ( |117| ). Here we 
have used the fact that a Bose distribution satisfies 



fix) 



1 

e x - 1 



= -[f(-x) + l] 



(119) 
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and energy and momentum conservation 

Pi + P2 = P3 + P4 



(120) 



However, while the fact that / is given by the locaLequilibrium Bose 



distribution in (117) ensures that C22H] — 0, one finds 
More precisely, the term in (|S9|) reduces to 



that C 12 [f] ? 0. 



(l + /l)/2/3-/l(l+/ 2 )(l + / 3 ) 



-/3[/I-Ai c -|ni(v„-v c ) 2 ] 



1 



(i + A)hk 



(121) 



using energy and momentum conservation. The expression in the square 
bracket in (121) only vanishes if the condensate and non-condensate are in 
diffusive equilibrium, which requires that the chemical potentials must be 
equal 

ft = fi c + ^m(v„ - v c ) 2 . (122) 

When we perturb the system, this may not be true, ie, the two components 
may be out of diffusive equilibrium. 

We can now derive hydrodynamic equations for the non-condensate by 
taking moment£_pf the Boltzmann equation, the standard procedure used in 
classical gases. The first moment gives a continuity equation with a source 
term: 

dn 



dp [Cf = C 12 [f)} - ^ = -V • (nv„) + r 12 [/], 



where we have used the fact that C22 [/] = and 

dp 



(123) 



(2tt) 3 
dp p 



(2tt) 3 m 
dp 



/(P,r,t) 
/(p,r,i) 



(2nf 



Ci 2 [/,$] 



ri 2 [/, 



More explicitly, we find [compare with ([71 
12[Jl ~ (2tt)5 [ 



(124) 



Canberra-rev: submitted to World Scientific on February 1, 2008 31 



J dpi j * dp 2 y dp 3 <5(mv c + pi - p 2 - p 3 ) 



x S(e c + £i-S2- e 3 )(l + h)hh 



T12 



-P\fi-Mc—im(v n 



1 



(125) 



We note that t\i is a collision time which describes the Cn collisions between 
the condensate and non-condensate atoms. Adding (123) to the continuity 
equation in (pT 



dn c 
dt 



-V-(n c v c )-r 12 [/, $], 



we see that the source term T 12 cancels out to give 
d(n c + n) 



dt 



= -V • (n c v c + nv n ). 



(126) 



(127) 



Thus our theory gives the correct continuity equation for the total local density 
n = n c + h. 

Similarly, one finds 



J dpp {Cf = C 12 [f}} - run + iw*) 

= -VP(r, t) - nVU(r, t) - m(v n - v c )r 12 [/], 



where the kinetic pressure is given by 



dp 

(2^y 



(p - m-v n ff(p,Y,t). 



(128) 
(129) 



Finally, the second moment gives 



j dpp 2 [Cf = C 12 [/]} - ^ + V • (Pv n ) 



.§PV.v„ 



17 



ri 2 [/]. 



(130) 



The detailed derivation of these equations is not important h ere. It involves 
straightforward manipulations using the explicit form of / in (117). 

The hydrodynamic equations ( |123| ), ( |128| ) and ( |130| ) describe the non- 
condensate in terms of three new "coarse-grained" local variables: 



h(r, i), v n (r, f) and P(r,t). 
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These are coupled to the two additional local variables which describe the 
condensate: 



n c (r,i), v c (r,t). 

We note that the two condensate equations of motion given by (^50|) are always 
"hydrodynamic" in form. In contrast, it is only in the collision-dominated 
region that the non-condensate dynamics can also be described in terms of a 
few collective variables. Both components exhibit coupled, coherent collective 
motions at the same frequency. This is the essence of two-fluid supcrfluid 
behaviour, familiar in liquid 4 Hc studies but still an unexplored frontier in 
trapped Bose gases. 

We will now discuss the linearized version of our condensate and non- 



condensate equations for local equilibrium as given by (cfj), (122), (128) and 



(130). We work to first order in the fluctuations around static equilibrium, 

n = ho + 5n, v n = <5v„, P = Pq + 8P 

n c = n c0 + Sn c , v c = <5v c . (131) 

What is new about the two-fluid hydrodynamic equations derived above is 
the role of the source term T^/]. In a linearized theory expanded around the 
static equilibrium Bose distribution /o (where r^f/o, $0] vanishes), one finds 

r i2 [/>] = STutf,*] = -fo%*5wt, (132) 

T 12 

where 

fidiff(r,t) = £(r,t) - M c (r,t). (133) 



Here r" 2 (r) is the C\ 2 collision time defined in ( |125| ) with both components be- 
ing in static equilibrium (see also (fT^)). The linearized coupled hydrodioiamic 
equations for the two components are given by the ZGN' equations LIS 

= -V • (n c0 5v c ) - ST 12 
m— -± = -V5fi c (134) 



and 

dSn 
~dt 



at 



= - V • (h 5v n ) + 5T 12 



mho—^- = - V5P - 5nVU (r) - 2gh a [VSn c + VSh] 
dSP 5 2 2 

"ST = - o V ' ( P 5 Vn) + T^ n • VP - ^Tltf ST 12 , (135) 

at 3 3 3 
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with ST12 given by (132) and (within the TF approximation) 



SfJ-c = g5n c + 2gSn. (136) 
We note that to lowest order, 5Ti2 does not appear in the second equation in 



( |135| ). In addition, in the last term in the third equation in (135), we have used 
McO — Uq = —gn C Q. The condensate couples into the non-condensate equations 
in ( |135[ ) directly via the mean field terms gSn c and, indirectly, through 5Ti2 
(see Section ||) . 

Despite appearances, this coupled set of hydrodynamic equations form a 



closed set. We have 9 scalar equations in (134) and (135), while they appear 
to involve 10 fluctuating variables 



5n c ,5v c ; 6n,5v n ; SP,Sfj,diff- 



(137) 



However one can showQ that 5/idiff can be written as a linear combination of 
the variables SP, Sn and Sn c and thus we are in fact left with only 9 variables. 
In Section 0, we prove that (134) and ( |l35| ) are equivalent to the two-fluid 
hydrodynamic equations of Landau when applied to a trapped Bose gas, but 
only when the condensate and non-condensate are in diffusive local equilib- 
rium, ie, jU c (r,t) = jl(r,t). It will turn out that the rate at which dfidiff(r,t) 
relaxes to zero is controlled by a relaxation .time t m related to the C12 colli- 
sions. The Landau two-fluid hydrodynamics 113 is only valid for low frequency 
oscillations which satisfy lot^ < 1, as we discuss in Section^. 



5 Two-fluid hydrodynamics in Bose liquids and gases 

To put our new two-fluid equations derived in Section ^ into some sort of 
context, we first briefly review the more general theory developed by Landau 
in 1941 to explain superfluidity in liquid 4 He. 

The original discovery of superfluidity in liquid 4 He is associated with the 
famous 1938 papers of Kapitza in Moscow and Allen and Misener at Cam- 
bridge. (I cannot resist remarking that Allen and Misener had been graduate 
students at the University of Toronto where they carried out (with Burton) 
the pioneering studies on vanishing viscosity in thejaeriod 1935-1937). These 
and subsequent experiments in the next few years c3 showed that superfluid 
4 He could exhibit very bizarre behaviour compared to ordinary liquids. This 
led to the development of a two-fluid theory of the hydrodynamic behaviour of 
liquid 4 He by Landau (1941). An earlier but less complete version of Landau's 
hydrodynamic equations was developed by Tisza in the period 1938-40. Tor 
further discussion of this early history, I refer to a recent article of mine. 13 
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In this early work, superfluidity (the term was coined in 1938 by Kapitza) 
was entirely associated with the relative motion of the_normal fluid and the 
superfluid components under a variety of conditions, ell The main point was 
that while the normal fluid exhibited finite viscosity and thermal conduc- 
tivity typical of ordinary fluid, the superfluid component (which exhibited 
irrotational flow) did not. In more recent times, the aspect of superfluidity 
which has been emphasized (see Ch.4 of Ref. l3 and Leggett's lectures in this 
book, for example) are those most directly tied to the fact that the superfluid 
velocity is associated with the gradient of the phase of the macroscopic wave- 
function <&(r, t). Given that this point of view is indeed more fundamental, it 
is still crucial to understand why superfluidity persists even in the presence 
of a dissipative normal fluid. This question can be addressed using two-fluid 
hydrodynamic equations. 

In essence, Landau developed his generic two-fluid hydrodynamics by gen- 
eralizing the standard theory of classical hydrodynamics to include the equa- 
tions of motion for a new "superfluid" degree of freedom. We recall that 
classical fluid dynamics was developed well before one knew about the exis- 
tence of atoms. Since the work of Maxwell and Boltzmann in the 1880's, we 
know the "coarse-grained" hydrodynamic description of a fluid in terms of a 
few quantities like n(r, t) and v(r, t) is only valid when-the collisions between 
atoms are strong enough to produce local equilibrium. Ea It only describes low 
frequency phenomena, where the condition in (^) is satisfied. 

In his 1941 paper, Landau did not connect the superfluid component 
with the motion of a ^'Bose condensate." Indeed, he rejected the efforts by 
Tisza and F. London B to use a Bose-condensed gas to get some insight into 
superfluid 4 He. However, since the period 1957-1965, Landau's superfluid 
degree of freedom has been understood microscopically in terms of the complex 
order parameter $(r, t). As noted earlier, the superfluid velocity field v s (r, t) 
is related to the gradient of the phase of $(r, t), as given by In the same 
period, it was also realized that there are two distinct kinds of quantum fluids: 

(a) Bose fluids (associated with a Bose condensate wavefunction $) 

(b) Normal Fermi fluids (associated with the key role of a Fermi surface). 

These t.wr» kinds of quantum fluids were magnificently described in two 
books LOO by Nozieres and Pines written around 1965, although the one 
on superfluid Bose liquids was only published in 1990 (it was in wide circula- 
tion as a preprint before then) . I think the clearest account of the connection 
between <&(r, t) and superfluidity in Bose fluids is still the discussion given in 
Chapters 4 and 5 of Vol II by Nozieres and Pines. Eil I highly recommend it 
to everyone in the BEC field. 
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Landau's pioneering work on superfluid 4 He has two separate aspects 
which are logically distinct but sometimes confused with each other: 

(a) The two-fluid equations describing hydrodynamic behaviour. These equa- 
tions are generic and apply (under certain conditions) to trapped Bose 
gases as well as to superfluid 4 He. In the period 1947-1950, Landau and 
Khalatnikov extended these equations to include hydrodynamic damping 
of the normal fluid (described by various kinds of viscosities, thermal con- 
ductivity, etc). This work is all described in the classic 1965 monograph 
by Khalatnikov. Il3 

(b) A "microscopic" theory of the elementary excitations describing the nor- 
mal fluid of liquid 4 He - the famous phonon - roton spectrum. Within his 
picture of a weakly interacting gas of quasiparticles, Landau and cowork- 
ers could calculate the thermodynamic and transport properties of su- 
perfluid 4 He. These quantities enter into the expressions given by the 
two-fluid equations for the first and second sound modes (velocity and 
damping). Of course, the roton part of the spectrum is not valid for a 
dilute Bose gas. 

As we have noted, Landau's original formulation of his two- fluid hydrody- 
namic equations was phenomenological in that the superfluid component was 
not given an explicit microscopic basis. The first derivation of the Landau 
hydrodynamic equations starting simply from the existence of the -macro- 
scopic order parameter <&(r,i) was given by Bogoliubov in 1963. E_ZJ This 
derivation built on Bogoliubov's earlier work on deriving hydrodynamic equa- 
tions for classical liquids without going through the intermediate stage of us- 
ing Boltzmann-like kinetic equations. While Bogoliubov's derivation is often 
viewed as being quite general, buried in his complex analysis is the assump- 
tion that the normal fluid and the superfluid are in local equilibrium with 
each other. Even today, probably the definitive account which formulates the 
various levels_Df theory for Bose superfluids is the classic paper by Hohenberg 
and Martin E3 published in 1965. It clearly shows (to the patient reader!) 
the central unifying role of <I>(r,i), summarizes the collisionless and hydrody- 
namic domains, and finally gives criteria for developing and judging various 
approximation schemes using Green's function techniques. 

With the preceding discussion as a preamble, I will now summarize the lin- 
earized two-fluid hydrodynamic equations of Landau, following the approach 
given in Ch. 7 of Ref. EJ. As noted, these equations are valid for both su- 
perfluid liquids as well as gases. The differences come in only at the last 
stage when we evaluate the thermodynamic coefficients appearing in these 
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equations, using the appropriate quasiparticle excitation spectrum. SEE 
The first two (linearized) Landau equations are familiar from ordinary fluid 
dynamics 

-fr + V • Sj = (138) 



= -V5P - 6nVV ex , (139) 

where we have included the effect of an external potential. Landau's work 
incorporated two components and hence the total mass density and mass 
current fluctuations are given by 

dp = mSn = 8p s + 5p n (140) 



mSj = Pso$v s +p nQ 5v n . (141) 

The new superfluid component was argued to only exhibit pure potential 
(irrotational) flow and carry no entropy. Thus Landau's final two hydrody- 
namic equations were new, 

m ^T = ~ vs " (142) 



^- + V • (s 6v n ) = .. (143) 
at 

Here the local equilibrium entropy density is s(r, t) = sq + 5s, the local equi- 
librium pressure is P(r, t) = P + 5P and the local chemical potential is 
n(r,i) = /-to + 6/j,(r,t). We note that the Landau hydrodynamic equations 
only involves 8 equations, in contrast to the 9 equations we derived at the 
end of Section ^. In particular, we see that Landau does not have separate 
continuity equations for the superfluid p s (r,t) and normal fluid p n (r,t) den- 
sities. Finally, it is assumed that these two components are always in local 
equilibrium with each other and hence the fluctuations are related by the 
thermodynamic identity 

n a Sp = -s Q ST + SP, (144) 
where no is the equilibrium total density. 
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For simplicity, we now specialize our analysis to a uniform super- 
fluid, where p s o, p„o, sq and other thermodynamic quantities are all position- 



independent. Multiplying (142) by n s0 gives [using (144)] 
d5v 



Pso- 



dt 



= -n s0 VSfi 



V(~s 5T + 6P) 

.^ V SP+^s V8T. 

no no 



Using this in (138) gives 
-V5P = 



Pso- 



d5v s 



<9<h 



dt 



PnO- 



dt 



n s0 
no 



VSP 



n s0 
n 



s VST + p n o 



dt 



(145) 



or 



PnO- 



d5v n 



dt 



n n p 
no 



VSP 



n s p 
no 



s VST. 



Combining ( |138| ) and ( |139| ) gives 
d 2 5p _ 



m V 



dt 



V 2 SP. 



Finally, (Q gives [using ( |146| )] 
d 2 Ss 



dt 2 



s V 



cMv„ 
~dt 



Po Po 
s d 2 Sp s 2 , 



PsO 
PnO 
PsO 
PnO 



V 2 (5T 



V 2 5T. 



(146) 



(147) 



(148) 

Po dt Po 

The last equation gives the local entropy fluctuations in terms of the local mass 
density bp and temperature ST fluctuations. Eq.( |148| ) can be re- written Ell in 
terms of the local entropy per unit mass s(r, t) = s(r, t) / p(r, t). Using 

—5s, (149) 
Po 



5s 



Po 



(148) takes on the simpler form 

d 2 5s _ 



PnO 
PsO 



V 2 <5T .. 



(150) 
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Expanding SP and ST in terms of Sp and Ss fluctuations 



dp 



as 



5s 



(151) 



we see that (|150|) and (147) reduce to two coupled scalar equations for Sp and 



Ss. Inserting the normal mode solutions 



Sp, Ss oc e 



i(k-r— ut) 



(152) 



one finds these coupled algebraic equations have two phonon solutions ui 2 
u 2 k 2 , where u 2 is the solution of the quadratic equation: 



dP 
~dp~ 



PsO 
PnO 



Ts 2 



T 

c v 

dP 
~dp~ 



1 dP 
~p~o~dT 

= 0. 



PsO 



Ts, 



PnO Cv 



(153) 



ds I 
dT I 



is the equilibrium specific heat per unit mass. 



Here c v 

The coefficients in ( |I53| ) are daunting, but only involve equilibrium ther- 
modynamic quantities which can be calculated for any given supcrfluid. We 
emphasize that ( |153| ) is valid for both a uniform Bose-condensed gas and 
superfluid 4 He, assuming that the superfluid and normal fluid are in local 
equilibrium with each other. However, as we shall show, the detailed char- 
acteristics and behaviour of the two phonon modes (firs±_and second second) 
are quite different in a Bose gas and in superfluid 4 He. cBeB 

A key feature about superfluid 4 He is that (typical of any liquid) 
dP/dT\ p ~ 0. This can be shown to be equivalent to C p ~ C v , where C PiV 
are t he specific heats at constant pressure or constant volume. In this case, 
(153) reduces to 



u 2 (A + B) + AB = 0, 



with two solutions u 2 = A, B : 

2 dP 
Ul = ^ 

2 PsO 

U 2 = 

PnO 



first sound 



Ts 2 



second sound. 



(154) 

(155) 
(156) 
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Working out the associated motions, one finds the first sound mode (oj = u%k) 
involves the in-phase motion of the superfLuid and normal fluid components. 
Moreover, one can show that first sound is essentially a pressure wave. In 
contrast, the second sound mode (u> = 112k) involves the out-of-phase motion 
of the two components, with 6j ~ or 

Pno5v n = -p s0 5v s .. (157) 

This corresponds to an almost pure temperature wave. The successful detec- 
tion in 1946 of a second sound mode (or teraapature wave) was of tremen- 
dous significance in low temperature physics. rWA The good agreement of the 
measured second sound velocity 112 with the Landau expression in ( |l5(j ), cal- 
culated using the postulated phonon-roton quasiparticle spectrum, vindicated 
both aspects of the Landau theory of supcrfluid 4 Hc. 

In contrast, in a (Bose-condensed) gas, the pressure and temperature 
fluctuations are strongly coupled and hence dP/dT\ p in ( |153| ) plays a signif- 



icant role. We recall that in a classical gas, one has C p /C v = |. Evaluating 
all the thermodynamic derivatives using our finite T microscopic model (the 
single-particle HF spectrum in (|5^) replaces the roton spectrum of superfluid 
4 He) and working to lowest order in gn c /kBT, one finds (after very lengthy 
calculations) 



2 



5 k B T g 5 /2(zg = 1) 
3 m #3/2(2:0 = 1) 

2 9 n c0 . 



0(gn) (158) 



(159) 



m 

It turns out that the first sound mode in a Bose gas is largely an oscillation 
of the thermal cloud (the normal fluid). The second sound mode is largely 
an oscillation of the condensate (the superfluid). Both involve density fluctu- 
ations and thus will have significant weight in the dynamic structure factor 
related to the density response function. The quite different features of first 
and second sound in Bose gases vs superfluid 4 He should be remembered when 



reading standard texts ahoujLsiinerfluid He. Results equivalent to (158) and 



(159) are given in Refs 



Comparing (159) with (|28|), it is clear that in a uniform Bose-condensed 
gas, the hydrodynamic second sound mode at finite T smoothly extrapolates 
to the Bogoliubov phonon mode in the collisionless region. The fact that 
the first sound velocity u\ does not depend on the interaction strength g to 
lowest order is typical of ordinary sound waves in any gas. However, one must 
remember that interactions (collisons) play a crucial indirect role in enforcing 
dynamic local equilibrium. 



Canberra-rev: submitted to World Scientific on February 1, 2008 40 



We now turn to a discussion of the generalized two-fluid hydrodynamic 
equations we derived at the end of Section || Again, for simplicity, we limit 
our analysis to a uniform Bose-condensed gas. If we can ignore vorticity in 
both fluids, we can introduce two velocity potentials 

<5v c = V0 c (r,i) 

6v n = V<f> n (r,t), (160) 



and reduce the equations in ( |134[) a nd (135) to three equations for tj> c , 4> n and 
S^diff [the latter is defined in ( |133| )]: 

d 2 4> c 9 2 a H 

m ~gt?~ = 9 n c(N 4>c + 2gn V <p n + —OfXdiff (161) 

m ^T= I !5 + 2<?n ) V 2 cp n + 2gn c0 V 2 4> c -l — Sfidiff (162) 
Ot z y3 n J 3 

^™ = \ 9 n c< N^ n - gn c0 V^ c - (163) 
or 3 t p 

Here r M is a new relaxation time governing how 5fi c uff relaxes to (ie, /i — > 
/i c ). It is found to be related to the C\i collision time r' defined in ( |8C| ) by 

- = —„ (164) 
whcj|c_the dimensionless hydrodynamic renormalization factor oh is given 



by 



J32 



_ 2 u 2J'"U (165) 

|P (1 - 7 ) + 2gh n c0 + §37o™c + ffl^o 

For completeness, we recall that 

Po(^o) = ^55/2(^0 = c-^-) 
^o(^o) = T3 33/2(^0) 

^° = Jrf-A** 1 ^' (166) 

where the Bose-Einstein functions are g n {z) = X^Si 2 an( l ^0 is the equi- 
librium thermal de Broglie wavelength. Calculation shows that cth becomes 
large as T — > Tbec an d thus r M can become very large near the superfluid 
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Figure 5. Squares of the first and second sound velocities (normalized by the first sound 
velocity of the ideal gas at T = Tg^c) vs - 1 /Tbec- The solid lines are for the Landau 
limit WTp <C 1 while the dashed lines are for o;r M 3> 1. The curves were generated for 
gn/kw^BEC = 0.2. The interesting crossover at low temperatures was first pointed out in 
Ref. Cj. It is reproduced by our equations, even thoughrjt-«ccurs when gn c Q ~ kgT which 
is outside the region of validity of our model. See Refs. LrCJ. 



transition. This is not in contradiction with the fact that the present dis- 
cussion is for the collision-dominated region [LOT22 <C 1, where T22 is some 
relaxation time associated with the C22 collision integral in (jS8[)1, 

It is clear from the fact that we have 3 coupled equations for the 3 variables 
4> c , <fi n and Sfidiff, we will obtain a new mode, in addition to the usual first and 
second sound oscillations discussed earlier. It is convenient to first eliminate 
S/^diff using the solution of (|l63| ), 



-iwSndiff = ofl« c o(-A: 2 )0„ - gn cQ (-k 2 )(fi c % - 
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or 



x . 2 2 

1 — (WT U O 



(167) 



Using this in (161) and ([163) gives 



mw C = gn c0 1 



1 — ILOT,, 



k <p c + 2gn 1 + 



Off "cO \ k 2 



3(1 - iujT^j n a 



5Po 
3 n 



2gh 



2oh n 



cO 



9(1 - iwT^) n§ 



fc 0„ 



2gn 



cO 



1 



(168) 



(169) 



3(1 - iwT^) n 

These two coupled equations for (p n and <f) c are easily solved and we obtain (as 
expected) first and second sound modes. Clearly the velocities U\ and u 2 will 
now depend on the value of cot^, although it turns out that this dependence 
is not very strong. This is shown by the results in Fig. [5]. 

The ojt^ — > limit is of special interest since the equations ( 168 ) and 



( |169| ) can then be shown to be completely equivalent to the predictions of the 
Landau two-fluid equations discussed earlier in this secti on. To be prec ise, 
the first and second sound velocities u\ 2 obt ained from ( 168 ) and Jl69|) in 



the limit wr M — > agree with those given by ( 153 ). This equivalence makes 



physical sen se si nce in the limit of r M — * 0, one sees that S/j,diff — * very 



rapidly [see (163) and (167)]. Thus we have proven that the Landau two- fluid 
hydrodynamic equations are correct if the condensate and thermal cloud are 
in diffusive local equilibrium. This, as noted earlier, is an (sometim es i mplicit) 
assumption in all previous derivations of the Landau equations. Lj'Lj 

This Landau limit {ojt^ — > 0) is, in fact, very subtle in the context of our 
microscopic calculation. We see that in this limit, there are still correction 
terms in ( |168| ) and ( |169| ) which are proportional to the hydrodynamic renor- 
malization factor an as defined in (165). These terms are crucial in ensuring 
that ( |168| ) and (169) reproduce the results of the Landau two- fluid equations. 
Another way of seeing this is that even though 5fj,diff — * when t m — > 0, we 
note that 6T12 in (132) is still finite. Using (167), one finds 

ST 12 = -n c0 a H - ^cpn) k 2 . (170) 

In conclusion, one might say that the hydrodynamic renormalization factor 
on contains a key part of the physics buried in the hydrodynamic two-fluid 
equations of Landau. 
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Figure 6. A schematic illustration of the predicted hydrodynamic mode spectrum for a uni- 
form Bose gas. (a) Above Tbec> on e has the usual first sound mode and a zero- frequency 
thermal diffusion mode, (b) Below TbeCi the standard Landau-Khalatnikov two-fluid hy- 
drodynamics predicts first and second sound, but no remnant of any zero-frequency mode, 
(c) In the ZGN' theory presented in Sections 4 and 5, one has ftjst and second sound modes 
as well as a zero-frequency relaxational mode. Based on Ref. Ft . 
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In the opposite limit ljt^ 3> 1 (which can arise near Tbec) we see all 



the terms proportional to <jh in ( |168| ) and (169) are negligible. This domain 
is missed in the Landau two-fluid equations. It describes situations in which 
the condensate and thermal cloud are out of diffusive equilibrium with each 
other. It would be of great interest to look for this kind of phenomenon in 
trapped Bose gases. 

One can work out the frequency of the new mode associated with the 
dynamics of SfAdiff and for a uniform gas, it is well approximated by 

Thus, in general, our two-fluid hydrodynamic equations predict the existence 
of a relaxational mode peaked at zero frequency. One can improve the the- 
ory to include deviations from local equilibrium [ie, deviations of / from / in 
( |116| )1. This involves a Chapman-Enskog kind of calculation familiar in the 
theory of classical gases and gives rise to hydrodynamic damping, tij The nor- 
mal fluid (non-condensate) equations of motion have new terms corresponding 
to shear viscosity (77) and thermal conductivity (k) transport coefficients E3e3, 
where we recall that k and rj are proportional to some T22 collision time and 
hence go as l/g 2 - The damping of first sound, second sound and the re- 
laxational mode due to small deviations from local equilibrium among the 
thermal atoms (C22[/, *&] )|-b|as been recently worked out in detail by Nikuni, 
Zaremba and the author. Ej Of particular interest is the effect on the relax- 



ational mode, which is now described by [compare with (171)] 

" 1 



h Auk 

T,. 



(172) 



Effectively the relaxation mode is strongly coupled into thermal conduction 
processes (or flow of heat). The expression for the coefficient A in (|l72| ) is 

somewhat complicated E3 but above Tbec (7 > 0), we find the mode reduces 

to the well-known thermal diffusion mode, 

ojr = -i-^r = -iD T k 2 .. (173) 
n Op 

This result strongly suggests that our new relaxational mode below Tbec is 
the renormalized version of the usual thermal diffusion mode above Tbec- 

We recall that in the standard two-fluid hydrodynamic theory of Landau, 
the thermal diffusion mode below Tbec disappears below TbeCi with, its 
spectral weight going into the emerging second sound doublet. In Fig. ^, we 
schematically illustrate the different hydrodynamic mode spectra predicted 
for above and below Tbec- 
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To illustrate the essential physics, the detailed discussion in this section 
has been limited to a uniform Bose-condensed gas. The analysis can, of 
course_hfi extended to trapped Bose gas. In particular, the out-of-phasc dipolc 
modeEju in the hydrodynamic limit is of special interest. This corresponds to 
centre-of-mass oscillation of the equilibrium condensate and non-condensate 
density profiles, 

n c (r,t) = n c0 (r -r) c {t)) 

n(r,t)=no(r-Tj n (t)), (174) 

with 

dVn dq c 
V "~ ;Vc = ^' (1?5) 

One finds that the hydrodynamic equations give (as expected) an undamped 
in-phase normal mode, with v„ = v c and lu — luq, where u>q is the trap 
frequency. This is a generalized version of the T = Kohn mode discussed 
at the end of Section ^[ In addition, however, there is an out-of-phase mode 
satisfying 

N c v c = -7Vv„, (176) 

with a frequency different from the trap frequency ujq. The frequency of this 
mode is shown in Fig. [?]. Calculation shows that this out-of-phase mode is 



only damped by its coupling to the relaxational mode given by (171). There 
is no hydrodynamic-type damping from the finite transport coefficients of 
the normal gas. Calculations are in reasonable agreement with the observed 
damping of such a out-of-phase mode. Ell A careful study of this out-of-phase 
dipole mode as a function of the temperature would be a nice way of probing 
the unusual hydrodynamics of a Bose-condensed gas. 



In these lectures, I havi 
experimental investigations 



e any attempt to analyze the available 
of collective modes at finite tempera- 
tures. These pioneering experimental studies are very promising but we need 
more systematic investigations, especially as a function of temperature and 
density. A key reason why such studies are perhaps a unique way of probing 
the many-body dynamics of superfluid gases c3 is that one can measure col- 
lective mode frequencies very accurately (at T = 0, with errors of only a few 
percent). 

In this section, we have put emphasis on deriving two-fluid hydrodynamic 
equations by starting from an approximate but still microscopic model. Such 
a derivation has been carried out recently by Zaremba, Nikuni and the au- 
thor for a trapped Bose gas. The pioneering work on deriving two-fluid 
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Figure 7. (a) Mode frequencies for the in-phase (solid dots) and out-of-phase (open dots) 
dipole modes vs. temperature for 2000 Rb atoms in an isotropic parabolic trap (see Ref. 13 
for values of the physical parameters used) . (b) Condensate (r] c ) and non-condensate (r) n ) 
amplitudes for the out-of-phase diptrk] mode, (c) Fraction of atoms in the condensate as a 
function of temperature. From Ref. E3. 



hydrodynamic .equations for a uniform dilute Bose gas was by Kirkpatrick 
and Dorfman. till In this regard, it is useful to point out that the linearized 
Landau two- fluid hydrodynamic equations [as given by (138) - (143)] are ex- 
pected to be exact as long as the two components are in local equilibrium. 
Thus, for a uniform system, the exact first and second sound velocities are 
given by the solutions of (153). The only question is how to calculate the 
various thermodynamic functions which are involved in this equation. The 
analogous Landau equations for a trapped gas allows us to go past the sim- 
plified microscopic models we have used in Section ^] to derive equations of 
this kind. In particular, one should be able to use the Landau two-fluid equa- 
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tions jas-a. direct probe of the superfluid density, just as one does in superfluid 
4 He. Eil'cH One might be able to detect small differences between the mag- 
nitude of the superfluid density n S Q and the condensate density n C Q (at the 
level of our model in Section ^, n S Q and n c $ are equal). In this connection, we 
note that in the formal zero temperature limit, the Landau two- fluid equa- 
tions reduce to two coupled equations for a pure superfluid (p s = p, p n = 
at T — 0). This limit has been used E-a to find a more accurate version of the 
T = quantum hydrodynamic equations discussed in Section ^. 

In conclusion, I hope I have given some insight into why the two-fluid 
hydrodynamics of trapped Bose gases has so much potential interest. I hope 
some of the young experimentalists attending this Summer School will take 
up the challenge to study this new frontier. In a sense, the natural next step 
after understanding the dynamics of a pure condensate (T = 0) is to study 
the two-fluid hydrodynamics of trapped gases (at finite T) since now one has 
two components which can execute coupled coherent collective motions. In 
contrast, the collisionless region at finite T as discussed in Section || does not 
seem as interesting since the thermal cloud has always such a low density (see 
Fig. ||). As a result, jineaji field effects produced by the non-condensate are 
never very important LJO, above or below Tbeci an d thus no new many-body 
dynamics emerges which is different than already found at T — 0. What is 
new is the different mechanisms of damping of collisionless condensate motion, 
as discussed in Sections || and [|. 
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